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Abstract 



LIGO observatories in Livingston, LA and Hanford, WA may detect gravitational waves 
emitted from coalescing binary systems composed of two compact objects. In order to 
detect compact binary coalescence (CBC) events, LIGO searches utilize matched filtering 
techniques. Matched filtering is the optimal detection strategy for stationary, Gaussian 
noise, however, LIGO noise is often non-stationary, non-Gaussian. Non-stationary noise 
result in an excess of false candidate events, commonly known as false alarms. This thesis 
develops the r^ test to reduce the false alarm rate for LIGO CBC searches. Results of the 
search for primordial black hole binary systems (where each object has less than IMq), in 
LIGO's Third Science Run (S3) is also presented. 

Results of the r^ test are shown for several LIGO CBC searches, including the binary 
neutron star searches in the Third and Fourth Science Runs (S3/S4), the S3/S4 primordial 
black hole searches, and the binary black hole search in the first three months of the Fifth 
Science Run (S5). The r^ test significantly reduces the false alarm rate in these searches, 
while only falsely dismissing a small fraction of simulated events. 



X 



1. Introduction 



At present time, a worldwide network of interferometric detectors are searching for grav- 
itational waves. These include LIGO (1) - which consists of two observatories with three 
detectors located in Hanford, WA USA (2km, 4km) and Livingston, LA (4km) USA; VIRGO 
(2) - a 4km interferometer located in Pisa, Italy; and GEO (3) - a 600m interferometer located 
in Hannover, Germany. 

This thesis contributes to the search of gravitational waves, specifically those emitted 
from coalescing binary systems composed of two compact objects (such as a neutron stars or 
black holes), by providing a new method the reduce the rate of false candidate events. The 
method is called the r^ test and will be introduced in chapter four. Results of the search 
for primordial black hole binary systems (where each object has less than IMq), in LIGO's 
Third Science Run (S3) is also presented in chapter five. This thesis specifically uses data 
from the LIGO interferometers, although methods described could be used with data from 
the other gravitational wave detectors mentioned. 

Chapter two introduces gravitational radiation, gravitational waves, astrophysical sources 
of gravitational waves, and the Laser Interferometer Gravitational Wave Observatory (LIGO). 

Chapter three describes data analysis for coalescing binary systems, including the the- 
ory of matched filtering, methods for detecting gravitational waves from coalescing binary 
systems, tuning searches for coalescing binary systems, and vetoes. 

In the fourth chapter, we describe methods to reduce false alarms in coalescing binary 
system searches, including the veto, and a test to further reduce false candidate events, 
known as the r^ test. The last section describes how the r^ test is implemented into the 
current searches of gravitational waves using LIGO data. 
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The fifth chapter summarizes the search for primordial black holes in the Third LIGO 
Science Run (S3), including the tuning of the search, results for the r^ test, and the result 
of the PBH binaries search in S3. 

The sixth chapter summarizes results for the r^ test in several other LIGO searches 
including the S3/S4 binary neutron star searches, S4 primordial black hole search, and the 
S5 binary black hole search (epoch 1). 

The analysis presented in this thesis is the fruit of my own work, in collaboration with the 
Compact Binary Coalescence (CBC) analysis group, which includes members of the LIGO 
Scientific Collaboration (LSC) and Virgo. The results are under review by the LSC, thus 
they are subject to possible revision before journal publication and they do not necessarily 
reflect the opinions of the LIGO Scientific Collaboration, the CBC group, or Louisiana State 
University. 



2. Gravitational Radiation 



This chapter describes gravitational waves, their sources, and an experiment presently 
conducted to detect these waves. In section 2.1, we will provide the basis of gravitational 
waves physics as described by Einstein's theory of General Relativity. Section 2.2 discusses 
the various sources of gravitational radiation, and in section 2.3 we describe detectors being 
operated today to discover gravitational waves directly, the Laser Interferometer Gravita- 
tional Wave Observatory (LIGO). 

2.1 Gravitational Waves 

Einstein's theory of General Relativity (4) can be summed up in this simple, yet powerful 
equation: 

G^, = SttGV, (2.1) 

known as Einstein's equation. It essentially tells us how the curvature of spacetime, coded 
in the Einstein tensor G^j^, reacts to the presence of matter /energy, coded in T^u, the energy 
momentum tensor, where G is Newton's constant (5). Since Gf^i, and T^j, are symmetric 4- 
tensors, this equation is actually a set of 10 coupled non-linear partial differential equations. 
Few exact solutions are analytically known due to the non-linearity of the equations. One way 
of gaining intuition about physical solutions to the field equations is to make approximations. 
For example, in the weak field limit we assume the gravitational field is very weak and the 
spacetime approximates that of Minkowski space. The metric, g^jy, can be written as: 

9nu = Vnu + h^u-, \hf,u\ « 1 (2.2) 

where rj^y is the flat Minkowski metric and h^^i, is the small metric perturbation, which will 
influence the geometry of space-time, and thus the motion of particles and how light travels. 
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The weak field limit itself transforms the non-linear Einstein equations in vacuum {Tfj,u = 0) 
into linear equations, with further simplifications using the transverse traceless (TT) gauge, 
therefore making Einstein's equation become a wave equation: 

For a a wave propagating in the z direction, h/^i, is: 

where h_|_(t) and hx(t) are the two independent polarizations of the gravitational wave. 
For gravitational waves being emitted from a binary composed of two neutron stars, the 
polarizations of the gravitational wave are given by (6): 

hJt) = — -^{TTGMff'^il + cos^ l) cos(2n ft - 2^0) (2.5) 
AG 

hy, (t) = - ^ixi-KGMff'^ cos i sin(27r/i - 20o) (2.6) 

where i is the inclination angle of the source, 0o is the initial orbital phase, / is the frequency 
of the gravitational wave, r is the distance from the detector to the binary system, M is the 
total mass of the binary system, M = mi + 1712, fJ, is the reduced mass of the system, fi = 
mim2 / (mi-|-m2). The frequency of the gravitational wave, /, evolves as the two neutron 
stars inspiral toward one another and is given by (6; 7): 

53/8^15/8 

^ " 87rG5/8/x3/8Mi/'i(t^ - t)3/8 (^.7) 

where c is the speed of light and tc is the coalescence time. 

Physically speaking, what is the effect of this gravitational wave passing by an object? A 
classic example is to imagine a ring of particles in the x-y plane. As the gravitational wave 
propagates along the z axis, h+ would stretch the distances between the particles along x 
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axis while shrinking the distances in the y axis and vice-versa. The hx polarization of the 
wave would have the same effect, although it would be rotated by 45°. This is shown in 
figure 2.1. 




Phase 7t/2 7t 3t^2 3k 

Figure 2.1: The effect of a gravitational wave on a ring of particles. The wave propagation 
is perpendicular to the ring. The top row shows the effect of the wave + polarized, while 
the bottom row shows the effect of the x polarization (8). 

What type of astrophysical phenomena can cause a gravitational wave? This is described 
in the following section. 

2.2 Astrophysical Sources of Gravitational Waves 

There are a large range of astrophysical phenomena that could produce gravitational waves 
strong enough plausibly detectable by ground based interferometers. An example of such 
detectors is the LIGO, where phenomena that produce gravitational waves within LIGO's 
sensitive band, the 40Hz-4kHz, could be detected in the near future if the signal is strong 
enough to be above noise and is distinguishable from false alarms. 

The first category of gravitational wave sources include those that emit brief, broadband 
signals, such as those produced from the collapse of a supernova or coalescence of binary 
systems, called hurst signals. In particular. Type II supernova collapses can yield strong 
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gravitational waves (h~10~^°). Another example of a burst source source includes those 
producing gamma ray bursts (GRBs), where sources are speculated to be from a variety of 
astrophysical phenomena ranging from black hole mergers to solar flares. 

A second category includes periodic signals, resulting from rotating stars, such as pulsars. 
Studies performed on the first discovered binary pulsar (PSR 1913+16) by Russell Hulse and 
Taylor in 1974 provide the best empirical evidence we have of gravitational waves indirectly 
to date (9). This was done by studying the orbital period of the binary system and accurately 
predicted the speeding up as predicted by general relativistic emission of gravitational waves. 

The third source of gravitational waves is a stochastic background, gravitational waves 
emitted from the beginnings of the early universe (cosmological stochastic background) or 
from a large an ensemble of unresolved astrophysical sources. 

Gravitational waves can also be emitted by compact binary star systems, composed of 
neutron stars or black holes. General Relativity predicts that the shrinking of the binary 
orbit of these systems is due to emission of gravitational radiation, with the gravitational 
waves being emitted at twice the orbital frequency of the binary. This process continues until 
the orbit begins to rapidly shrink (inspiral) and gravitational waveform becomes a "chirp" 
signal, increasing in frequency and amplitude. To give a sense of the time scale for this 
process, imagine an inspiralling system composed of two 1.4 M© traversing through LIGO's 
sensitive band with a starting frequency of 40Hz, the signal would spend about 25 seconds 
in LIGO's sensitive band. In the case of low mass binary systems, we can approximate the 
waveform using the 2"*^ order Post-Newtonian expansion (10; 11). We will return to the 
discussion of coalescing binary systems and methods used to detect these signals in chapter 
3. 
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2.3 Laser Interferometer Gravitational Wave Observa- 
tory (LIGO) 

The Laser Interferometer Gravitational Wave Observatory (LIGO) is a joint Caltech and MIT 
project supported by the National Science Foundation. The observatory is composed of two 
sites, one in Hanford, WA and the other in Livingston, LA. The LIGO Hanford Observatory 
(LHO) contains both a 4km (HI) and 2km (H2) interferometers which are independent of 
one another, yet share the same vacuum system. The LIGO Livingston Observatory (LLO) 
contains one 4km interferometer (LI). Each of the interferometer arms are perpendicular to 
one another. The orientation of the Hanford Observatory is closely aligned to the Livingston 
Observatory, as to give a common response to a gravitational wave signal, the North-South 
arms of both detectors are along a great circle. This is depicted in figure 2.2. 




Figure 2.2: The two LIGO sites (14). 

Beyond the two LIGO sites, an international team of scientists has been created in order 
to analyze the data and to contribute to the improving of the experiment. This group is 
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called the LIGO Scientific Collaboration (LSC) (15). It includes scientists working at the 
LIGO laboratory sites and scientists representing over 52 institutions from around the world, 
including Louisiana State University. 



2.3.1 Design 

The LIGO detector design is a power-recycled Michelson interferometer with Fabry-Perot 
arm cavities, which have a set of equal length perpendicular arms of 4 km. The laser light 
source (Nd:YAG laser, operating at a wavelength of 1.06 fim with a maximum output power 
of low) is sent to the beam splitter which directs the light along the two perpendicular 
arms, which first pass through a partially transmitting input mirror (ITM) which forms a 
Fabry-Perot cavity with the end mirrors (ETM). The light from both arms recombine and 
interfere with itself destructively at the beam splitter, with the non eliminated light hitting 
the photodetector (AS port), this in fact is the signal we will use to detect gravitational 
waves. A simple version of the LIGO detector is shown in figure 2.3. 

ETMY 



ITMY 





RM 




laser 









ITMX ETMX 
photodetector 



Figure 2.3: A simplified schematic of LIGO interferometer. RM is the recycling mirror, 
ITMX,Y are the intermediate test mass (partially transmitting mirrors), ETMX,Y are the 
end test mass (mirrors). 
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The data taken which measures the strain signal while interferometer is in "lock" is used 
to determine if a gravitational wave has been observed. To acquire "lock" means to take 
the interferometer from a state where the mirrors move freely to a state where the light is 
resonant in the arm cavities and the power-recycling cavity, and the interference at the AS 
port is kept at the dark point. A more detailed explanation of the LIGO interferometers can 
be found in (16) and how to "lock" the interferometer in (17). The data taken in the locked 
state and when the feedback systems keep the detector in a "low noise" mode is referred to 
as "science data". This strain signal, s(t), is defined as: 

sit) = ' (2.8) 

where Ly.^y are the arm lengths on the interferometer. This signal is limited by the amount 
of noise in the interferometer, n{t), which will be described in the following section. 



2.3.2 Noise Sources 

There are several types of noise sources which could affect LIGO and limit LIGO sensitivity. 
These noise sources contribute to the measured output, denoted as h{f). We show in figure 
2.4, the three main noise sources which determine LIGO's design sensitivity. These are 
seismic noise, thermal noise, and laser shot noise. Seismic and thermal noise affect the 
interferometer by altering the differential change in the arm lengths measured in the anti- 
symmetric port, and thus the strain signal defined in equation 2.8. 

Seismic noise limits the strain sensitivity at frequencies f 40 Hz. There are several 
constituents that contribute to the measured seismic noise, including earthquakes (0.03 ^ f 
<^ 0.1 Hz), ocean waves (0.1 ^ f ^ 0.35 Hz), and noise made by the everyday living of people 
(1.0 Hz ^ f 3 Hz). This noise is reduced by combination of passive and active isolation 
techniques (18). 
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Thermal noise may be the limit in the 40 Hz ^ f ^ 200 Hz regime and is due to thermal 
fluctuations of the individual elements that compose the test masses (fused silica), including 
the coating on the mirrors, the substrates, and the steel wires that hold the optics in place. 

Laser shot noise affects the strain signal sensed in equation 2.8 and dominates for f ^ 200 
Hz. This is due to fluctuations in the stored power in the interferometer arm Fabry-Perot 
cavities from fluctuations in the vacuum field, alternatively interpreted as photon counting 
statistics (shot noise) on the light. 

The noise sources described above limit the overall sensitivity of LIGO. Actual LIGO 
detectors are limited by these and a variety of other noise sources, but the overall measured 
noise currently similar to the LIGO design goal. This is illustrated in figure 2.4. 




Frequency (Hz) 



Figure 2.4: Noise limits for the LIGO interferometer design (20). 



A more detailed description of all noise sources in the LIGO interferometers can be found 
in (19). 
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2.3.3 Science Runs 

As described in section 2.3.1, science data is taken when the interferometer is in a locked 
state and is in a low noise mode. Having both LIGO sites operating in a locked state results 
in times where coincident data can be taken. Once the LIGO sites achieved a sensitivity 
that was better than previous bar detectors, coincident data taking began, which alternated 
with commissioning to improve the sensitivity . This was organized into long periods of time 
of coincident data taking operation, known as science data taking runs, or science runs. The 
LSC has analyzed the data for gravitational waves in the data taken from these science runs 
offline. Up to the writing of this thesis, there have been 4 completed science runs with a 
current science run S5: 

1. SI August 23, 2002 to September 09, 2002 

2. S2 February 14 2003 to 14 April 14 2003 

3. S3 October 31 2003 to January 09 2004 

4. S4 February 22, 2005 to March 24, 2005 

5. S5 Began November 4, 2005 and is currently in progress. 

Through S4, no direct evidence of gravitational waves have been detected (21; 22; 23; 
24; 25; 26; 27; 28; 29; 30; 31; 32; 33; 34; 35; 36; 37; 38; 39; 40; 41; 42; 43). In S5, the 
experiment reached the goal in sensitivity where we are now surveying hundreds of galaxies. 
The improvement in strain sensitivity over the course of the science runs is shown in figure 
2.5. 
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10^ 10^ 10^ lo'* 

Frequency (Hz) 



Figure 2.5: Best strain sensitivities for the LIGO 4km Interferometers (45). 



3. Data Analysis for Coalescing Binary 

Systems 



In this chapter I introduce the fundamentals of gravitational wave data analysis for searching 
coalescing binary systems. We begin with the theory of matched filtering, followed by a 
description of the method used in a multi-detector search by the LSC in a multi-detector 
search, the "inspiral pipehne" . The chapter also discusses methods to tune inspiral searches, 
and vetoes used, both instrumental and signal-based. 



In searching for gravitational waves with well understood waveforms, such as "chirp" wave- 
forms from inspiralling neutron stars, matched filtering (46; 47) is the optimal detection 
strategy in Gaussian noise. The measured detector's stain amplitude is given by 



where n{t) is the detector noise (assumed in this chapter to be Gaussian distributed) and h{t) 
is the gravitational wave signal. An interferometric detector, such as LIGO (section 2.3), is 
sensitive to a linear combination of the two gravitational wave polarizations (equations 2.5, 
2.6), where the gravitational wave signal has the form: 



3.1 Optimal Matched Filtering 



s{t) = n{t) + h{t), 



(3.1) 



h{t) = F+h+{t) + F^h^{t) 



(3.2) 



with 



(1 cos^ 6') cos 20 



(3.3) 



X — 



COS 9 sin 2(f) 



(3.4) 
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where F+ and Fx are the two antenna response functions of the detector, 6*, and are the 
spherical coordinates of the sources sky position with respect to axes defined by LIGO's 
arms. 

Wc transform the template, h{t) (equations 2.5, 2.6, 3.2, 3.3, 3.4), into the frequency 
domain using the the second order post- Newtonian stationary phase approximation{48) 

hs{f)^ihc{f) (3-6) 

where / is the frequency of the gravitational wave, M is the total mass of the binary system, 
M — (mi + m2), ^JL is the reduced mass of the system, = mim2/(mi + 7712), and r] is the 
ratio ji/M. The template h{f) (units of 1/Hz) is a hnear superposition of hc{f) and hs{f), 
defined in equations 3.5, 3.6. The phase evolution of the chirp signal, ^'(/; M, 77), is given by 
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^'(/;M,77) = 27r/t,-20o-7r/4 + — + -r;ja;-3-167ra; 



/3715 55 



12877 

/ 15293365 , 27145 , 3085 2\ -1 



+ + 77 + 77^ x-^ , (3.7) 

V 508032 504 ' 72 ' ' ' ^ ' 



„/, -V3 



. = f^) (3.8) 

where 0o is the coalescence phase, and tc is the coalescence time of the binary. A combination 
of equations 3.5, 3.6, 3.7, and 3.5 yields the form of the chirp signal used in the search. An 
example chirp signal for a 1.4—1.4 Mq binary with a starting frequency of 40Hz is given in 
figure 3.1. 

Optimal performance is achieved by filtering the data s{t) in the frequency domain 
weighted inversely by the power spectral density of the noise. The complex output of the 
Wiener optimal filter is given by: 

<')-^r^'^j§l^'"'''df (3.9) 
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Figure 3.1: An example chirp waveform for a 1.4—1.4 Mq coalescing binary system, / = 
40Hz at t = Os, / = 1532Hz at t = 25s. A zoom of the first 0.5 seconds is shown in the 
bottom panel. 
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where s(/) is the Fourier transformed data of the detector and Sn{f) is the power spectral 
density of the detector noise. The goal is to locate the maxima of the output of the matched 
filter, 11-2(^)11, over the arrival time and phase, comparing the maximum value with the 
expected distribution of values in the absence of the signal. Furthermore, the output of the 
optimal filter for each template has a characteristic amplitude: 

we can then normalize the filter output, defining a signal-to-noise ratio (p): 

P(t) = ^ (3.11) 

A template has a time series while a candidate event (or trigger) is related to a local 
maximum of and is assigned an SNR: 

(3.12) 

(7 

Given a candidate with a signal-to-noise ratio (SNR) p, we can also infer another physical 
characteristic of the system, the effective distance: 

Deff = (IMpc)^, (3.13) 
P 

which is the distance at which an optimally oriented binary inspiral would produce the 
observed signal-to-noise ratio. Des is related to the physical distance of the binary, r, by 
the orientation of the binary with respect to the plane of the detector and the detectors two 
antenna response functions by (5): 

Deff = , ^ = (3.14) 

^Fl{l + cos2 i)2/4 + Fl C0S2 L 

A true inspiral signal would have a narrow peak in p(t) at the coalescence time, tc, the end 
of the template used, or the time at which the two stars begin to merge. If no gravitational 
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wave is present, all we would have is stationary detector noise and {p^if)) = 2, since p{t) is 
a random quantity, distributed with 2 degrees of freedom. 

Instrumental glitches in the data can lead to large p, even in the absence of a signal: 
these are called "false alarms". Further tests are needed to differentiate real signals from 
false alarms; these include coincidence between detectors, the veto, and the r^ test, which 
will be further investigated in 4. 

A more detailed explanation of how matched filtering is implemented in LIGO CBC 
searches can be found in (6). 

3.2 Searching for Gravitational Waves from Coalescing 
Binary Systems: The Inspiral Pipeline 

The search for gravitational waves from binary systems in a multi-detector search done by the 
LSC uses several steps, collectively called the "inspiral pipeline". The pipeline begins with 
the raw data from a set of interferometers and ends with a number of candidate events, which 
have passed a multitude of tests. These include signal thresholds, multi-tier coincidence, 
signal based vetoes, and data quality cuts. The pipeline itself can be run with simulated 
signals, either physically injected into the interferometer (hardware injections) or into the 
data itself (software injections), and provide a great test of the efficiency to detect inspiral 
waveforms. 

The inspiral search pipeline can be divided into the following stages. 

1. Data Collection: For the LIGO searches, we use a channel named DARM_ERR. This 
channel is sampled at 16384 Hz and is the error signal of the feedback loop, controUing 
Lx - Ly, keeping the differential arm length (equation 2.8) constant. 

2. Template Banks: The template bank generation begins by reading in the DARM_ERR 
data, down-sampling it to 4096 Hz, and applying a high-pass filter. The data is di- 
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vided into 2048 second "blocks". These "blocks" are further divided into 256 second 
segments, which are then overlapped with each other by 128 seconds, illustrated in 
figure 3.3. This is done to insure no corrupted output is included into the analysis 
due to edge effects of the Fast-Fourier- Transform. The average power spectral density 
of the noise, <S'(/), is also calculated, as well as the instrumental response function, 
which is used to calibrate the data and the power spectrum. Once this is complete, the 
template waveforms are generated for each "block" of data. Note that we require more 
than one template to scan to full parameter space, therefore using a bank of templates. 
The number of templates in a bank is governed by allowed mismatch in the signal with 
the template waveform, which ensures we lose < 3% in SNR of the signal. In the S4 
searches (39), the template banks per block of data was 3500 for BNS, 4500 for PBH, 
and 1200 for BBH. 

3. Matched Filtering I: As described previously in section 3.1, the templates generated 
from step 2 are run through each block of data. The SNR time series, is 
calculated for each template in each segment, where a SNR threshold, p*, is applied. 
Inspiral "candidates" are kept if p > p*.The SNR threshold (p*) is chosen low enough 
to allow a background of false alarms. 

4. First Coincidence: At this stage, the inspiral triggers from multiple interferometers are 
combined, and we look for triggers coincident in time and mass within some accuracy. 
The triggers that pass this test are kept, while the others are discarded. Demanding 
coincidence between interferometers significantly lowers false alarms. 

5. Triggered Template Bank: A new template bank is generated using only the parameters 
of triggers which survived the previous step. Their mass parameters are used to create 
the new bank. 
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6. Matched Filtering II: The templates generated from stage 5 are run through the data 
again. 

7. Vetoes: Several instrumental (sec. 3. 4) and signal-based (sec. 4.1, sec. 4. 2) vetoes are 
applied. 

8. Second Coincidence: We look again for coincidences between candidates left after in 
step 7. These are our candidate events. 

The pipehne is run on a given data set using Directed Acychc Graph (DAG), which 
describes the workflow, and is executed using the Condor High Throughput Computing 
System (49), which manages the execution of tasks. The DAG reads in a configuration 
file which contains locations of executables, science data lists, and search parameters. The 
output of each stage of the pipeline after the data collection are in the extensible markup 
language (XML) format, which stores in a tabular format the parameters and information 
about the triggers. The software used to perform all steps of the analysis and to construct 
the DAG is in the package LALAPPS{50). 

We also repeat the pipeline with time shifts to find the background rate of the false 
alarms. Time shifts are done by taking one set of interferometer triggers at one detector 
and shifting them in time, usually on the order of seconds. These time shifted triggers 
are then run through the first coincidence stage where any coincidence between the time 
shifted triggers and real interferometer triggers are now false alarms. The pipeline is also 
repeated for software injections (signals simulated in software) to determine the efficiency of 
the search. A visual representation of the pipeline described can be found in figure 3.2. 

3.3 Tuning the Search for InspiraUing Binaries 

Searches for gravitational waves from coalescing binary systems must be optimized using 
variety of methods in order to increase the likelihood of detecting gravitational waves. In 
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Figure 3.2: Illustration of the Inspiral Coincidence Analysis Pipeline (44). 
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Science segment Time 




1 28 s overlap 

Figure 3.3: Illustration of how single interferometer data is divided (22). 

other words, parameters used in different stages of the search pipeline must be tuned. Data 
from the playground, which constitutes roughly 10% of all the data from a science run taken 
in 600 second intervals every 6370 seconds, is used to do this. The tuning presented here will 
be primarily focused on what was used for LIGO's 3^'*^ and 4*'^ science runs for the binary 
neutron star (BNS) and primordial black hole (PBH) searches. 

A first step in tuning our searches is allowing a time coincidence window taking into 
account the precision in recovering the gravitational waveform's inferred coalescence time, 
tc, and the light travel time between detector sites. The precision is measured by injecting 
simulated gravitational waves into the archived data and running matched filtering. Since 
we know what the parameters of the injected waveforms are, we can take the output of the 
matched filters to see how well the detected parameters agree with the injected parameters. 
We histogram the difference between injected and recovered tc- The width of the histogram 
gives the timing accuracy. The S3 BNS search used a window of 4 milliseconds to find an 
injection in a single detector, for coincidence between sites we also include the light travel 
time, which is on the order of 10 milliseconds. 

Other coincidence windows are tuned. For non-spinning waveforms, the masses are the 



22 



main parameters, and we can therefore choose a coincidence test based upon the chirp mass. 
The chirp mass is given by: 

and the symmetric mass ratio of the binary is: 

Precision in the recovered chirp masses depends on the variance in the detector noise (equa- 
tion 3.10) and discreteness of the template bank. For the S3 and S4 PBH and BNS searches, 
we chose a window for both 77 and M.c based upon histograms the injected and detected 
parameters. 

A third test can be only done in the case of the two co-located IFO's, such as HI and H2 
by tuning the precision in the amplitude, or the effective distance, defined in equation 3.13, 
for candidates found in both IFO's. The precision in effective distance is given as 

2|HlDefr - H2Defr| 

where K is the value tuned with the software injections. In S3, the sensitivity of of HI and 
H2 differed significantly, and the test used was: 

iHlDeff - H2Defrl 

H1P,„ < . (3.18) 

The amplitude consistency check is not done between the Hanford and Livingston detectors 
since the detectors are not perfectly ahgned due to the curvature of the Earth; therefore we 
could not expect a similar amphtude between sites. 

3.4 Vetoes 



So far I have described the detection strategy used for detecting inspiral signals, the analysis 
pipeline, and how to tune such a search. I now bring focus to the characterization of vetoes. 
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Conceptually, a veto is a time where by instrumental inspection or some statistical calcula- 
tion, the data is atypical. There are three classes of vetoes; [1] instrumental vetoes - caused 
by problems found in the interferometer, [2] statistical vetoes - defined by correlations be- 
tween the gravitational wave channel and other channels, and [3] signal based vetoes, which 
determine if candidate events are consistent with the waveform. 

3.4.1 Instrumental Vetoes 

An instrumental veto occurs at a time when the instrument was known to malfunction or is 
sensing an environmental excitation and the malfunction is known to generate false alarms. 
In the ideal case, once the times due to all instrumental problems are vetoed, the distribution 
of remaining false alarms is as predicted for Gaussian data. 

The goal in veto investigations is several-fold: it should provide a high efficiency of 
dismissing false signals (high veto efficiency), while keeping the amount of science mode 
time used at a minimum (low deadtime), and it should have a large percentage of veto 
triggers that eliminate one inspiral event or more. It is a balance between these factors that 
we try to achieve, with the more analyzable time we have with Gaussian noise data, the 
greater the chance of detection. 

There are times during any of science runs when environmental disturbances and non- 
optimal interferometer performance times should be excluded from the analysis, even when 
the interferometer is in science mode(51). Large transients in the data causes most, if not 
all of the templates in the bank to ring off, and thus create false alarms, even in coincidence. 
Members of the LSC doing detector characterization develop data quahty (DQ) flags. Some 
of the tuning is done with candidates from the first matched filtering stage of the analysis 
pipeline (sec. 3. 2). These triggers are used as a primer to go further into the IFO performance 
and look for causes of these occurrences. The result of such work are data quality (DQ) flags, 
used to let the scientists know what was discovered in the IFO while it was in science mode. 
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A DQ flag can be created for a number of reasons, ranging from earthquakes affecting tlie 
instrument, dips of the hght intensity in the arms, airplanes passing by, and an unrelated 
channel triggering. It should be noted that not all DQ flags are created from inspiral triggers, 
most are in fact created from known problems with the interferometer. Some DQ flags do not 
produce false alarms in compact binary coalescence (CBC) searches, thus, tuning is required. 

3.4.2 Statistical Vetoes 

Problems in the instrument can be identified in some statistical ways . An example includes 
looking for correlations between "kleineWelle triggers" (52) on the gravitational wave (GW) 
channel and other auxiliary channels. kleineWelle is a wavelet analysis algorithm. An exam- 
ple of such a veto is the "AS_1 veto" (53), where times when there are kleineWelle triggers 
in both AS_I and AS_Q, with a ratio that makes the transient a very unlikely gravitational 
wave. 

3.4.3 Signal Based Vetoes 

Finally, there is another category of vetoes which includes signal-based vetoes, which are 
used to discriminate a real signal from background noise in a given search. An example of 
this is the veto and the r^ veto, which is based upon the time series of x^. The following 
chapter will describe these both in detail. 



4. Methods to Reduce False Alarms in 
Coalescing Binary System Searches 



This chapter introduces a veto currently implemented in searches of gravitational waves from 
coalescing binary systems, the veto, and introduces a new veto, the r^ test. In the last 
section, I will describe how the test is incorporated into the inspiral pipeline used in LSC 
searches. The results presented in this chapter are my own while collaborating with the CBC 
group. 

4.1 The Veto 

As described in Sec.3.1, matched filtering provides the optimal detection strategy for de- 
tecting inspiral waveforms in Gaussian noise. Noise artifacts present in non-Gaussian noise 
data generate a signal-to-noise ratio (p) large enough to exceed the given search threshold 
(p*) resulting in a false alarm. In order to further reduce these transient signals from being 
misidentified as real signals, we check whether the waveform has the expected accumulation 
of p in several time-frequency bins, illustrated in figure 4.1. This is analogous to breaking one 
template into several sub templates where each sub template produces p/p of its parent 
template's SNR p. 

For each candidate, we define the time series as (12): 

X\t) = ^ E \ziit) - Z{t)/p\' =pj:\p^{t) - p{t)/p\' (4.1) 
^ 1=1 i=l 

where p is the number of bins (sub templates), pit) is the matched filter (SNR time 
series) for the signal, and pi{t) is the SNR time series of the sub templates, i. If the noise 
in the detector was Gaussian noise, we expect x^(^) to be x^ distributed with 2p-2 degrees 
of freedom in the absence of true signals. The reduced x^ test threshold (x*^) for Gaussian 
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Time 

Figure 4.1: A visual representation of statistic. The solid black line represents the chirp 
waveform, while the gray boxes illustrate each time-frequency bin, p = 8 in this example 
(13). 



noise is: 



2 2 

*2 , X _ X 



(4.2) 



where the quantity of the right has an expected value of 1. In LIGO inspiral searches, we 
re-normalize to have an expected value ~ 2. 



r 
p 



(4.3) 



This quantity is referred to as r^ (7). For inspiral searches with LIGO data, a discrete 
template bank is used that causes a potential mismatch between the template waveform and 
signal (~ 5% or less). The mismatch causes the to scale with p^. The r^ must be modified 
to account for this. Therefore, we add a non-central parameter Sp^ to r^: 

= (4-4) 

where 5 is the mismatch between the template and the signal. This is called the weighted 
X^, or In order allow room for these mismatched signals, we apply a threshold on 
denoted ^*^. We therefore ask when p exceeds p*, that the triggers also have both: 



X^ 
p + Sp'' 



*2 



(4.5) 
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Since this is a new constraint on the outputted triggers, the test now becomes a veto, 
hence the veto. 

An example comparing the output for a high p simulated inspiral injection's and 
weighted time series is shown in figure 4.2. Notice in the zoomed in bottom panel in 
figure 4.2 where the threshold for the loud injection would have vetoed the calculated r^. 

4.2 A Test to Further Reduce False Alarms 

In the inspiral pipeline (section 3.2), triggers are kept during the second matched filtering 
stage if they have a SNR above a threshold (p*) and the weighted below a threshold (^*^). 
In Gaussian noise, the probability of these triggers being a false alarm can be estimated. In 
real data collected by LIGO for example, this estimate is seriously fiawed due to the presence 
of excess noise contributed from instrumental transients. This is exemplified by comparing 
segment of data's filter output for a single template's pit) and weighted x^ time series in 
simulated Gaussian noise and actual detector noise (LHO 4km data). This is plotted in 
figure 4.3. In the simulated Gaussian noise case, both time series stay relatively fiat with 
some small fiuctuation, while in the detector noise, we see times when the SNR rises and 
falls and the weighted x^ time series changing by a significantly larger amount. Figure 4.4 
shows histograms of the SNR and the weighted x^ time series. Notice here the SNR time 
series overlaps in both cases, with the LIGO data having a tail. The weighted x^ has a slight 
shift in overlap, due to the presence of instrumental glitches. We want to develop a method 
to characterize the excess noise seen in these figures, and therefore create a veto to reduce 
the background. This can be done based of the time history of SNR p{t) or weighted x^ 
time series, or some variant thereof. 

Take as an example figure 4.5. The time series, p(t), of the false alarm and simulated 
inspiral injections behave differently up to the coalescence time. The weighted x^ time series 
has greater fiuctuations for the false alarms up to the inferred coalescence time, making it 



28 



Simulated Inspiral Injection 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



inferred coalescence time (s) 



10^ 



10^ 

CD 
CO 



10^' 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 
inferred coalescence time (s) 



10^ 





CO 
> 

10^ 



-0.1 -0.05 0.05 0.1 0.15 

inferred coalescence time (s) 



Figure 4.2: A simulated inspiral injection SNR and weighted time series using LIGO 
data. The top panel is a plot of the SNR time series, the middle panel plots the weighted 
time series and x^/P- The bottom panel is a zoom of the middle panel. The dashed lines 
represent the set thresholds. 
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Figure 4.3: The SNR and weighted time series for a 128 second segment of LIGO data 
and simulated stationary Gaussian noise 

easier to distinguish it from a true signaL Therefore, we choose to use create a test based 
upon the time series as a method to search for excess noise. We use r^ (equation 4.3) 
instead of the weighted (equation 4.4), due to the smooth behavior of r^ and also since the 
weighted has a minor effect until we reach the coalescence time of the simulated inspiral 
injection as shown in the bottom panel of figure 4.2. The test will measure the consistency of 
the candidate signal using r^(t) by scanning a given amount of time prior to the coalescence 
time at a given threshold. This test is otherwise known as the r^ test. 

The r^ test will require two parameters, the r^ threshold, r*^, and a window of time prior 
to the coalescence time to search. At,,,, which is on the order of seconds. The r^ test looks 
up to the time before the inferred coalescence time of the trigger since the templates do not 
include the merger waveforms. The test will determine how much time the r^ time series 
stays above r*^ in At*, with the result denoted as At. A visual representation of the test 
is given in figure 4.6. An example of how the test performs for an ensemble of simulated 
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Figure 4.4: A histogram of figure 4.3. The top panel is a histogram of the SNR time series. 
The bottom panel is a histogram of the weighted time series. 
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Figure 4.5: A comparison of SNR and weighted time series for a false alarm (top two 
panels) using LIGO data and a simulated inspiral injection (bottom two panels) in simulated 
stationary Gaussian noise. The dashed lines represent the set thresholds. 
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Figure 4.6: Illustration of the r^ test. 
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inspiral injections and time slide triggers (false alarms) for the S4 binary neutron star search 
are shown in figure 4.7. 

In figure 4.7 we see two distributions. A distribution of simulated inspiral waveforms 
that begin to separate from the background time slide triggers (false alarms) at SNR = 10, 
and a distribution of false alarms that separate from the simulated waveforms for SNR < 10 
and At > 0.0002 seconds. The contours plotted in the figure denote two separate methods 
to veto these false alarms while falsely dismissing few simulated inspiral injections. The first 
is a dashed black contour which includes a bottom left rectangular region that is defined by 
two parameters. The first is defined as a segment of constant time duration called r^ (in 
the case shown in the figure, r^ = 0.0002) which begins at the SNR threshold of the search 
and extends to the vertical dashed black line, which is the second parameter. The second 
parameter is the segment of constant SNR extending to form the long side of the rectangle 
called pr (in the figure, pr = 10). Using this black dashed contour line as a veto for this 
particular data set where any points within the rectangular regime are eliminated results in 
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Figure 4.7: test example result from the S4 binary neutron star search from background 
triggers (time slides) and simulated inspiral injections. We chose to plot At plus 10 millisec- 
onds to allow false alarms and simulated inspiral injections having At ~0 visible in the plot. 
The dashed black line has the chosen parameters r^ = 0.0002 and pr = 10. The magenta 
line has the chosen parameters of rc = 0.007 and rp = 1.1 . 
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35% of the false alarms being vetoed while 0.001% of the simulated inspiral injections are 
falsely dismissed. The use of the r^ test's results as a veto is called the r^ veto. The black 
dashed contour described in this example is called the veto vl. 

The second is a magenta dashed contour includes rj, Pr, and two additional parameters. 
The two additional parameters describe the power law region of the contour: 

/(p) = re X p^- (4.6) 

where rc and rp are parameters chosen to create the slope. Using this magenta dashed 
contour line (with ic = 0.007, rp = 1.1), where any points above will be eliminated for this 
particular data set results in roughly 36% of the background time slide triggers being vetoed 
while falsely ehminating 0.001% of the simulated inspiral injections for the figure given. The 
magenta dashed contour described in this example is called the veto v2. 

Note the four parameters can be tuned according to a given data set and inherently so 
the fraction of vetoed false alarms and falsely dismissed triggers will change. This will be 
shown in chapters 5 and 6. The next section describes how this test could be implemented 
into an inspiral pipeline. 

4.3 Implementation of the Test Into the Inspiral 
Pipeline 

The r^ test is implemented into the inspiral pipehne whenever the is calculated. The 

r^ test is incorporated into the inspiral search pipeline (Sec. 3.2) in the second matched 
filtering stage. The r^ test pipeline is described as: 

1. Characterize the candidate trigger by the signal to noise ratio (p), keep candidates 
with p > p*. 
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2. Measure the consistency of the candidate signal with the expected chirp waveform, 
keeping candidates with x^/(p + Sp^) < ^*^. 

3. Compute the time series x^/p, count the time the trigger stays above a given time 
threshold, r^, in a given interval. At*, prior to the inferred coalescence time, tc, of the 
trigger. Store the calculated value as At. 

The last step is what has been added to the inspiral pipehne. The following chapters will 
include examples of results from implementing the pipeline described. 



5. The Search for Primordial Black Holes in 
LIGO's 3rd Science Run 

A black hole composed of mass < 1.0 Mq is believed to be a primordial black hole (PBH) 
(29), since it could not have been created as a product of stellar evolution. These compact 
objects may have formed in the early, highly compressed stages of the universe immediately 
following the big bang. PBH's are also theorized to be candidate Massive Astrophysical 
Compact Halo Objects (MACHOs), which could be some fraction of the dark matter in the 
Large Magellanic Cloud. A binary system composed of two PBH's will emit gravitational 
waves that may be detectable by LIGO (5). This chapter describes the search for primordial 
black hole binary systems in the third science run, including the tuning of several coincidence 
parameters, how the r^ test was employed, and the result of the PBH binaries search from 
the S3 run. The results are included in (39). The results presented in this chapter are the 
author's while collaborating with the CBC group. 

5.1 The Third LIGO Science Run 

The third LIGO science run was conducted from October 3, 2003 to January 09, 2004. All 
three LIGO detectors at the two observatories were in operation. The best sensitivity curve 
for each IFO is shown in figure 5.1. 

Of the numerous binary systems that could be searched for gravitational wave emission, 
primordial black hole (PBH) binary systems were chosen to search for this thesis. The masses 
for the PBH search ranged from 0.35 Mq to 1.0 Mq. We used a low frequency cutoff, fi, of 
lOOHz. The longest waveform duration, L'max, was 22.1 seconds while the average number 
of templates used per block of data was 4500. The range of masses for the PBH search were 
from 0.35 M© to 1.0 Mq. An example PBH template waveform is given in figure 5.2. 
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Figure 5.1: LIGO S3 Run best sensitivities for all interferometers (54). 
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Figure 5.2: A sample PBH waveform, / = lOOHz at t = Os, / = 2023Hz at t = 22s. A zoom 
of the first 0.5 seconds is shown in the bottom panel. 
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In order to search for these signals, we use a data analysis pipeline as described in section 
3.2. The first step is to determine to amount of data that will be used by collecting the times 
when the three interferometers were in science mode. We also incorporate separate segment 
lists which include times where the quality of the data is non-optimal for each of the detectors. 
Once this information is collected, we calculate the times when two or more IFO's are in 
operation, this will represent the analyzed times to be searched and is summarized in table 
5.1. Since second coincidence is essential to the search, a viable candidate event must have at 
least two detectors observing the same event. We required both of the Hanford detectors to 
be in operation, due to the possibility of scattered hght from one of the out of lock detectors 
interfering with the other detector. 

Table 5.1: Times analyzed when at least two detectors were operating. The times in paren- 
theses exclude playground times. 

S3 

H1-H2-L1 times 184 (167) hrs 
H1-H2 times 604 (548) hrs 
Total times 788 (715) hrs 

5.2 Tuning the Search 

In order to tune a search for inspiralling binary systems (as described in section 3.3), several 
parameters must be tuned. This is done is order to correctly characterize and optimize the 
analysis pipeline for detection. Recall that two thresholds are introduced into the pipeline 
at various stages in order to reduce the number of noise transients, these are the SNR and 
the test thresholds. We chose p* = 6.5, = 10.0, with 16 bins. Recall p{t) and 
are calculated using equations 3.11 and 4.5. 

For the PBH binary system search, the x^ test (12) provides a measure of the quality-of- 
fit of the signal to the template. We can also combine the SNR (equation 3.11) and the x^ 
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value of the event, defined as the effective SNR, p^g, by 




where p is the number of bins used in the test; the specific value of the parameter 
250 is chosen empirically. This parameter was chosen in order to effectively separate our 
simulated inspiral injections with the background triggers (time slides) shown in figure 5.4, 
since using the SNR of each coincident inspiral injections or background triggers did not 
effectively separate the two as shown in figure 5.3. We expect pcs ~ p for real signals with 
relatively low SNR, and low effective SNR for noise transients with high x^- We can also 
assign to each candidate in coincidence a combined SNR, pc, defined by 

N 

(Pc)|BH = EPeff,i (5.2) 
i 

where i is the detector index. 



• H1L1 False Alarms 

• HI H2 False Alarms 




Figure 5.3: SNR of coincident time slides (false alarms) and coincident simulated inspiral 
injections for the S3 PBH search. 
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Figure 5.4: The effective SNR statistic for the S3 PBH search. 



Several coincidence parameters are also tuned. These include the coincident time between 
triggers (AT), the coincidence in the chirp mass (AAic), and the coincidence in t] (called 
Ar/). These parameters were tuned using results from injections. The results are plotted in 
figure 5.5 for AT, figure 5.6 for AJUc, and figure 5.7 for Arj. The chosen parameters based 
on these figures are given in table 5.2. 

Table 5.2: Summary of the S3 PBH coincidence windows. The second column gives the 
coincident-time windows column; we need to take into account for time of flight between 
detectors (10 ms between LI and H1/H2 detectors). 



AT (milliseconds) 


AMc (Mq) 


Ar] 


4x2 


0.002 X 2 


0.06 



As described in section 3.3, we can tune the error in the effective distance (equation 
3.18) for simulations done on the H1-H2 coincident triggers by comparing the background 
time-slide triggers as shown in figure 5.8, a value of k = 0.45 was chosen. 
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Figure 5.5: Timing difference (AT) histogram. 
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Figure 5.6: Mc difference (AAic) histogram. 
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Figure 5.7: 7] difference [At]) histogram. 

Finally, once all the various parameters have been tuned, we run the analysis pipeline 
from start to end and cluster the final coincident triggers. They are clustered by locating the 
trigger with the largest SNR within a 22 second window, which is the length of the longest 
template. 

5.3 Results for the Test 

The r^ test as described in section 4.2 was computed in the S3 PBH search. The parameters 
chosen for the test are given in table 5.3. 

The results of the r^ test for the S3 PBH search are plotted in figure 5.9. We see here 
a clear distinction between the time slides background triggers and simulated PBH binary 
system injections. The solid black line was chosen as a parameter region to veto any points 
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Figure 5.8: The H1H2 effective distance cut for the S3 PBH search. The bottom panel is a 
histogram of the H1H2 effective distance cut. The sohd black line is the value of k chosen 
{k = 0.45). 
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Table 5.3: Test Parameters for the S3 PBH search. 



r*^ At^, (seconds) pr I'd (seconds) Tc 
15 2.0 13 2 X 10"^ - - 




Figure 5.9: Resuh of the r^ test for the S3 PBH search using the parameters given in table 
5.3. The solid black line is the r^ veto chosen, where points within it being eliminated. 
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within. Using the sohd sohd black hne as a veto therefore results in eliminating 26.5% of 
the time shde background triggers and while falsely dismissing 0.001% of the injections. 

The result proves the power of using such a test on a binary inspiral search, where 
previously, coincidence between triggers and the test were the strongest vetoes available. 
And most importantly, it brings the possibility of finding gravitational waves in LIGO I and 
future detectors that much closer. In chapter 6, we will see the r^ test apphed to several 
other science searches and how the veto parameters were chosen. 

5.4 S3 PBH Search Result 

A cumulative histogram of the combined statistic, (Pc)pbh ^ defined in equation 5.2, of the 
loudest coincident triggers in the S3 PBH search is shown in figure 5.10. We estimate the 
background in the following way: 

1. We shift the time of the LI detector by 5 seconds from the true time, and the time of 
the H2 detector by 10 seconds from the true time. 

2. We look for coincident triggers, which are false coincidences (thus false alarms) due 
to the shifted time. 

3. We use the same pipehne (section 3.2) used for the zero time shifted analysis, including 
all vetoes. 

4. The resulting distribution of coincident triggers is considered the background of false 
alarms. 

5. We repeat the procedure for 50 different time shifts, and construct a cumulative 
histogram of the mean number of false alarms versus combined statistic, and also 
calculate the standard deviation in each bin. 
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No triple coincident foreground candidate events or background events were found. We 
show in figure 5.10 the foreground, or double coincidences found, and the expected back- 
ground. The number of double coincidences found was consistent with the measured back- 
ground. 

Cumulative histogram of Number of events vs Statistic 
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n 
E 



'" 70 80 90 100 110 120 130 

Combined Statistic 

Figure 5.10: Result of the S3 PBH search. The cumulative histogram of the combined 
statistic, (pc)pBH) coincident candidates events (A) and the estimated background from 
time slides (+) with la deviation ranges for the S3 PBH binaries shown (shaded region). 
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6. The Test in Other LIGO Searches 



In this chapter we revisit the r^ test, as described in chapter 4, apphed to several other LIGO 
searches for gravitational waves. This includes the search for binary neutron stars (BNS) in 
S3/S4, the search for primordial black holes (PBH) in S4, and the search for binary black 
holes (BBH) in the first three months of S5, also known as epoch 1. Each section will include 
a brief summary of the individual search. The searches in this chapter were done in the CBC 
group, where the author was a critical member on the discussions of the r^ test application. 

6.1 S3/S4 Binary Neutron Star Searches 

The LIGO's third science run was conducted from October 31, 2003 to January 09, 2004 
and the fourth science run ran from Feb 22, 2005 to March 24, 2005, with all all three LIGO 
detectors at the two observatories in operation. 

The first implementation of the r^ test was done on the S3 binary neutron star search. 
The BNS searches (39) used post-Newtonian templates in the range 1.0 Mq < mi,m2 < 
3.OM0, which for the S3/S4 searches are summarized in table 6.1, where iVb is the number 
of templates in the bank and -Dmax is the length of the longest template. The coincidence 
parameters are listed in table 6.2. 



Table 6.1: Target Sources of the S3/S4 BNS Search. 



^min 



max 



(Mo) /i(Hz) iV, Ana.(s) 



53 BNS 

54 BNS 



1.0 
1.0 



3.0 70 2000 10.0 

3.0 40 3500 44.4 



The r^ test parameters in the S3/S4 BNS search are given in table 6.3. The results of the 
test are given in figures 6.1 and 6.2, where the black line denotes the r^ veto chosen. The 
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Table 6.2: Summary of the S3/S4 BNS Coincidence Windows. 



AT (milliseconds) 


AMc (Mo) 


Ar] 


S3/S4 BNS 5x2 


0.01 X 2 


0.10 



use of the veto significantly lowers the background, vetoing about 43% of the false alarms 
for S3 and 35% for S4, even after the veto and other vetoes. The veto falsely dismisses 
only 0.001% of the simulated inspiral injections for both runs. 

Table 6.3: r^ Test Parameters for the S3/S4 BNS Search. 

r*^ At^, (seconds) r^ (seconds) r^ rp 

S3/S4 BNS 10 2.0 10 2 x 10"^ - ~ 




Figure 6.1: Result of the r^ test for the S3 BNS search using the parameters given in table 
6.3. The solid black line is the r^ veto chosen, where points within it being eliminated. 
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Figure 6.2: Result of the test for the S4 BNS search using the parameters given in table 
6.3. The solid black line is the r^ veto chosen, where points within it being eliminated. 



Table 6.4: r^ Test Results for S3/S4 BNS Searches 



Falsely Dismissed Injections (S3,S4) Vetoed False Alarms (S3,S4) 



0.001, 0.001 



43.0, 35.0 



6.2 S4 Primordial Black Hole Search 
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The S4 PBH searches (39) use post-Newtonian templates in the range 0.35 Mq < mi, 1712 < 
I.OMq, which is summarized in table 6.5. The coincidence parameters chosen are included 
in table 6.6. 



Table 6.5: Target Sources of the S4 PBH search. 



mmin(M0) 


"^max(M0) 


/l(Hz) 






0.35 


1.0 


100 


4500 


22.1 



Table 6.6: Summary of the S4 PBH coincidence windows. 



AT (milliseconds) 


AM, (Mq) 


Ar] 


4x2 


0.002 X 2 


0.06 



The r^ test was done on the S4 primordial black hole search (39) . The parameters are 
given in table 6.7. The results of the test are given in figure 6.3, where using the black hne is 
the r^ veto chosen and results in gives a significant lowering of the background on the order 
of 35% after the veto and other vetoes, while falsely dismissing 0.001% of the simulated 
inspiral injections. 

Table 6.7: r^ Test Parameters for the S4 PBH Search. 

r*" At,^ (seconds) f)^- (socoiids) r^. ip 
~S4 15 2.0 13 2 X 10"^ - ~ 



6.3 S5 Binary Black Hole Search (Epoch 1) 

S5 started on November 4, 2005. We analyze the "first epoch", from November 4, 2005 
to February 6, 2006, with all all three LIGO detectors at the two observatories were in 
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Table 6.8: test Results for S4 PBH Search 



Falsely Dismissed Injections Vetoed False Alarms 
% 0.001 3K0 
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operation. This search targeted black hole binaries of component mass in the range 3.0 
M© < mi, 777.2 < 3O.OM0, with of total mass no higher than 35.0 M©. For the S5 epoch 
1 search, we decided to use standard post-Newtonian templates with a test. A sample 
waveform is shown in figure 6.4. BBH simulated inspiral injections were done from 4 different 
waveform families which included Pade-Tl, Taylor-Tl, Taylor-T3, and Effective-One-Body 
(55). 

Since the test was used, we could therefore incorporate the r^ test as well. The 
parameters for the r^ test is given in table 6.11. The results of the test are given in figure 
6.5. Note this was the first time we decided to use a non-zero value of rp. In figure 6.5 we 
see an overlap between the false alarm triggers and simulated inspiral injections, possibly 
due to the shorter length of the BBH templates. This results in not being able to veto most 
false alarms beyond a given SNR as we saw in section 6.1, but does significantly lower the 
background on the order of 19%. 

Table 6.9: The target sources of the S5 BBH search. 

(^o) /l(Hz) iVb Anax(s) 

3.0 25.0 40 300 7.0 



Table 6.10: Summary of the S5 BBH coincidence windows. 

AT (milliseconds) /\Mc (M©) 
20 X 2 5.0 X 2 



Table 6.11: r^ Test Parameters for the S5 BBH search. 



r*^ At* (seconds) r^ (seconds) r^ rp 



2 X 10-^ 0.025 1.05 



S5 10 



6.0 9 
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Figure 6.4: A sample BBH waveform, / = 40Hz at t = Os, / = 701Hz at t = 7s. A zoom of 
the first 0.5 seconds is shown in the bottom paneL 



Table 6.12: r^ test Results for S5 BBH epoch 1 



Falsely Dismissed Injections Vetoed False Alarms 



% 



0.12 



19.1 
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SNR 



Figure 6.5: Result of the test for the S5 BBH search using the values given in table 6.11. 
The solid black line is the r^ veto chosen, where points within it being eliminated. 



7. Conclusions 



We have devised a new test, that we call the test, that greatly reduces the number of 
false alarms in gravitational wave searches of coalescing binary systems. We have shown the 
actual efficiency of this test using data from the LIGO interferometers from several science 
runs. The reduction in the rate of false alarms ranged for example from 43% for the search for 
binary neutron star systems in the third science run, to 19% in the search for the coalescence 
of binary black hole systems in the first three months of the fifth science run, while achieving 
a low false dismissal rate of simulations ranging from 0.001% to 0.12%, respectively. The 
test will be incorporated into future LIGO searches such as the searches for binary systems 
in the fifth science run. A search for primordial black hole binary systems (where each object 
has less than IM0) in LIGO's Third Science Run (S3) was conducted where results from the 
number of double coincidences found was consistent with the measured background. 
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